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ABSTRACT 

We present a kinematical study of the optical ejecta of GK Per. It is based on proper motions 
measurements of 282 knots from ~20 images spanning 25 years. Doppler-shifts are also computed 
for 217 knots. The combination of proper motions and radial velocities allows a unique 3-D view 
of the ejecta to be obtained. The main results are: (1) the outflow is a thick shell in which 
knots expand with a significant range of velocities, mostly between 600 and 1000 km s""'^; (2) 
kinematical ages indicate that knots have suffered only a modest deceleration since their ejection 
a century ago; (3) no evidence for anisotropy in the expansion rate is found; (4) velocity vectors 
are generally aligned along the radial direction but a symmetric pattern of non-radial velocities 
is also observed at specific directions; (5) the total iJa+[N II] fiux has been linearly decreasing 
at a rate of 2.6% per year in the last decade. The Eastern nebular side is fading at a slower 
rate than the Western one. Some of the knots displayed a rapid change of brightness during the 
2004-2011 period. Over a longer timescale, a progressive circularization and homogenization of 
the nebula is taking place; (6) a kinematic distance of 400±30 pc is determined. These results 
raise some problems to the previous interpretations of the evolution of GK Per. In particular, 
the idea of a strong interaction of the outflow with the surrounding medium in the Southwest 
quadrant is not supported by our data. 



Subject headings: ISM 
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jets and outflows; stars: individual (GK Persei) — novae, cataclysmic variables; 



1. INTRODUCTION 

Nova Persei 1901 (GK Per) was discovered by 
Anderson in 1901 February 21 at a visual magni- 
tude of 2.7. Only one day before, its brightness 
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must have been fainter than 12 mag ([Williams 
19011 ). It reached the maximum light on Febru- 



ary 22 when it was as bright as Vega (0.2 mag). 
GK Per was the first object around whic h super- 
luminal light echoes were detected (e.g. iRitchev 



19021 ). The actual ejec ta of the nova o utburst was 



discovered in 1916 by Barnard ( 1916[ ) and is still 
easily observable now: indeed, it is considered to 
be the longest lived and most energetic nova rem- 
nant ever found. 

The earliest optical images of the remnant 
showed an emission only i n the Southwest (SW ) 



side of the remnant (e.g. Seaauist et al.l [l989l) . 
Gradually more material became visible around 
the central source, however, the asymmetric na- 
ture of the remnant has remained. Radio observa- 
tions obtained during 1984-1986 showed that the 
SW quadrant of the remnant is the source of non- 
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therm al (synchrotron) emission (jSeaauist et al 
This would be caused by particle accel- 
eration and enhanced magnetic fields in a fast 
shock where the expelled matter is meeting some 
pre-existing material, most probably the result of 
previous mass ejections from t he system, possibl y 
an ancient planetary nebula ( Bode et al. l2004[ ). 
This scenario is supported by the analysis of the 
X-ray emission from the remnant (Balman 2005). 
A lumpy, highly asymmetric X-ray nebula is ob- 
served, with brightness enhancement in the SW 
quadrant and wing-like extensions in the South- 
east (SE) and Northwest (NW) directions. It is 
concluded that a density gradient in the circum- 
stellar medium causes the expansion of the X-ray 
shell to be faster in the NW and SE directions, 
reaching velocities as high as ~2600 - 2800 km s~^, 
than toward the SW, where velocities are less 



than half of that, around ^1100 km s 



This 



strong interaction is expected to cause decelera- 
tion of the expans ion of the SW part of the outflow 
(iBode et al.ll2004l ). 



At a distance of half a kiloparsec (e.g. iMcLaughlin 
1963) and with an expansion velocity of a thou- 
sand km s^^, the apparent growth of the nebula is 
about one arcsec per year, which is easily resolv- 
able from ground-based optical imagery even on a 
timescale of few months. Therefore GK Per offers 
the opportunity to study the evolution of a stellar 
outflow in real-time, and to determine not only 
the expansion velocities but also their derivatives 
(i.e. acceleration, thus the forces at work), or de- 
viations from pure radial expansion, an important 
information that is generally not available for dy- 
namical studies of astrophysical outflows. If the 
expansion in the plane of the sky is coupled with 
Doppler shift velocities, information about flve of 
the six dimensions of the phase-space is gained, 
and an almost complete dynamical picture can be 
drawn. For these reasons, starting on 2004 we 
embarked in a program of frequent imaging mon- 
itoring of the expansion of the GK Per remnant. 
We present here the results of the work. Doppler 
shift velocities of 217 knots from long-slit spec- 
troscopy are combined with high-precision mea- 
surements of the motion in the plane of the sky 
for 282 knots. The latter ones are computed us- 
ing good-quality, ground-based ifa-|-[NII] images 
spanning 25 years, with a particularly intensive 
time coverage since 2004. 



At the beginning of the article we describe our 
data and archival data. Then designations used 
in this article are given (Section \^ . Proper mo- 
tions and radial velocities for individual knots are 
discussed in Section [S] and HI and the kinemati- 
cal analysis is presented in Section [T] Section |8] 
presents the brightness evolution of the nebula, 
and Section |9] presents the discussion. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

2.1. New Imaging 

Most of the imaging was carried out at the 2.5m 
Isaac Newton Telescope (INT) at the Observatorio 
del Roque de los Muchachos (ORM), La Palma, 
Spain. Images were taken with the Wide Field 
Camera (WFC) and a narrow band Ha fllter hav- 
ing a central wavelength of 6568 A and a bandpass 
of 95 A. The pixel scale is 0".33 pix^i. 

Other images were taken with the Andalucia 
Faint Object Spectrograph and Camera (AL- 
FOSC) at the 2.6m Nordic Optical Telescope 
(NOT) of the ORM. An Ha fllter (6577/180 A) 
was used, and the pixel scale was 0".189 pix^^. 
Both the INT and NOT fihers are broad enough 
to include the emission from the [Nil] 6548 A and 
6583 A doublet, which arc the strongest optical 
emission lines in the ejecta, and to cover its whole 
radial velocity range (±1000 km s~^). 

All images were obtained with a seeing bet- 
ter than 1.3 arcsec full-width at half-maximum 
(FWHM). 

2.1.1. Archival imaging 

We also searched for archival image s of the neb- 
ula. INT images obtained on 1999 (iBode. et al 



2004 ) with the WFC, and on 1987 (Scaauis t et al 



19891 ) with the Prime Focus Cone Unit (PFCU, 
pixel scale 0".74 pix~^) were downloaded from the 
Isaac Newton Group archive. The 1999 data have 
been obtained with the same fllter as our recent 
INT-hWFC observations. The 1987 frame was ob- 
served in a broadband R filter. 

The [Nil] Hubble Space Te lescope (HST) 1995 
and 1997 images discussed by Shara et al.l (2012) 
were downloaded from the MAST public archive. 
They were obtained with the Wide Field Planetary 
Camera 2 (WFPC2): the relatively large remnant 
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fills the field of view of the camera, and is sampled 
at a 0".l pix-i by the WF detectors, and 0".045 
pix~^ by the PC CCD. In both epochs, the filter 
F658N was used, which has a central wavelength 
6590 A and a bandpass of 22 A. This is narrower 
than the range covered by the ejecta and some 
knots with the most extreme radial velocity shifts 
are missed. 

The Digitized Sky Survey (DSS) images from 
the Palomar 48-inch Samuel Oschin telescope 
(Pal48) obtained on 1953 (pixel scale 1".7 pix'^) 
and 1989 (pixel scale 1".0 pix~^) were retrieved 
via the Aladin interface. In both cases obser- 
vations were done on red-sensitive photographic 
plates with a bandpass approximately equivalent 
to an R filter. In 1953 a plate Kodak 103a-E and 
in 1989 a IIIaF was used. The 1953 image is used 
for visual comparison with the recent images, but 
not used for any quantitative measurements. The 
1989 image could instead be used to measure the 
position of several knots recognized in our most 
recent imagery. 

The log of all the imaging observations is pre- 
sented in Table [1] 

2.1.2. Imaging data reduction 

The basic reduction (bias, flat field, combina- 
tion of sub-exposures, and cosmics removal) of 
most of the images was done using standard rou- 
tines in HST data were reduced using 
tasks in the package dither. All the frames were 
cropped to the same field of view (FOV) of 6.4x6.4 
arcmin^ around the central star, except for the 
data from 1987, 1995, and 1997 which originally 
had a smaller FOV. 

The next step was to carefully register all im- 
ages to the same astrometric reference frame us- 
ing stars from USNO-Bl.O catalogue with proper 
motion errors smaller than 0".002 yr~^. About 
sixty stars were used in every frame, except in the 
older frames which had worse seeing and smaller 
FOV. Our first INT image, obtained on 20040129, 
was chosen to be the reference frame. The Python 
based Kapteyn Package was used to get the ref- 
erence frame FITS header world coordinate sys- 



IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 



tem applied pixel-by-pixel to the reference frame. 
The other images were pixel-by-pixel matched to 
the reference frame using tasks wcsxymatch, ge- 
omap, and geotran in IRAF. The latter put all the 
frames in the same pixel scale of the reference im- 
age, 0".33 pix^^. Results were checked using some 
twenty stars in the proximity or inside the rem- 
nant. The standard deviation of the position of 
these stars from the reference frame in the X and 
Y directions is 0.08 arcsec except for the images 
taken before 2004 where it is occasionally slightly 
larger. We adopt this figure as representative of 
the uncertainty introduced by the astrometric reg- 
istration of images of different epochs. 

To investigate brightness variations of the rem- 
nant, all the INT-I-WFC frames marked with * in 
Table [U were also flux matched. First, all the 
frames were degraded down to seeing of the ref- 
erence frame, 1.2 arcsec. Secondly, frames were 
flux matched using forty-eight well exposed stars 
in the FOV. No stars inside the remnant were used 
to avoid possible contamination from the moving 
knots. The rms of fluxes of the reference stars on 
different frames is 1.3 %. 

2.2. Spectroscopy 

Medium-resolution, long-slit spectra of the neb- 
ula were obtained during two observing periods. 
In 2007 January, the INT-hlDS (Intermediate Dis- 
persion Spectrograph) with grating R1200Y and 
a 1".5 slit was used. This configuration covers 
a spectral range from 5730 A to 7610 A at a 
dispersion of 0.47 A pix . In 2007 September, 
the NOT-I-ALFOSC with grism #17 and 0".5 slit 
was used, giving spectral range from 6350 A to 
6850 A, and a dispersion of 0.26 A pix~^. The 
spectrographs slit was oriented at a dozen of dif- 
ferent position angles with the aim to cover as 
many knots as possible. Details of the spectro- 
scopic observations, including the orientation of 
the slit, the exposure times, and the spectral reso- 
lution with the adopted slit widths, can be found 
in Table [21 These data were reduced and wave- 
length calibrated using the standard routines in 
the package longslit in IRAF. 

3. DESCRIPTION OF THE NEBULA 

The optical nebula of GK Per, displayed in 
Figure [1] has a roundish, knotty morphology. 
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Table 1 
Log of the imaging data 



Date 


JD 


Tclcscopc+ 


Filter 


Total cxp. 


Number 


Original seeing 






Instrument 


CW/FWHM (A) 


time (sec) 


of frames 


(arcsec) 


1953-01-10 


2434387.62 


Pal48 


phot, plate 


2400 


1 


5.0 


1987-07-27" 


2447003.69 


INT-I-PFCU 


R 


300 


1 


2.8 


1989-10-05 


2447804.87 


Pal48 


phot, plate 


5100 


1 


4.9 


1995-11-08*' 


2450030.42 


HST-fWFPC2 


6590/22 


1400 


2 


0.2 


1997-01-08*' 


2450457.18 


HST-I-WFPC2 


6590/22 


1200 


2 


0.2 


1999-11-30*" 


2451513.56 


INT-fWFC 


6568/95 


7026 


3 


2.3 


2004-01-29* 


2453034.37 


INT-fWFC 


6568/95 


900 


3 


1.2 


2004-08-04* 


2453222.66 


INT-fWFC 


6568/95 


1500 


5 


1.2 


2004-09-03* 


2453252.67 


INT-I-WFC 


6568/95 


1500 


5 


1.1 


2004-12-30* 


2453370.43 


INT-I-WFC 


6568/95 


1200 


4 


1.1 


2005-09-16* 


2453630.72 


INT-(-WFC 


6568/95 


900 


3 


1.0 


2006-10-12* 


2454021.53 


INT-I-WFC 


6568/95 


1200 


4 


1.0 


2007-01-04* 


2454105.49 


INT-I-WFC 


6568/95 


300 


1 


1.0 


2007-03-01* 


2454161.42 


INT-fWFC 


6568/95 


2100 


7 


1.2 


2007-09-03 


2454347.68 


not+alfosc 


6577/180 


1800 


3 


0.6 


2007-09-05 


2454349.71 


NOT-I-ALFOSC 


6577/180 


1800 


3 


0.5 


2007-09-07* 


2454351.74 


INT-fWFC 


6568/95 


1200 


4 


1.1 


2008-01-08* 


2454474.42 


INT-fWFC 


6568/95 


1200 


4 


0.8 


2008-08-19* 


2454698.69 


INT-fWFC 


6568/95 


4650 


31 


1.1 


2008-09-12 


2454722.70 


INT-fWFC 


6568/95 


1200 


4 


1.3 


2008-11-14* 


2454785.49 


INT-I-WFC 


6568/95 


1500 


5 


1.2 


2009-02-11 


2454874.43 


NOT-I-ALFOSC 


6577/180 


1200 


2 


0.8 


2009-12-08 


2455174.48 


INT-I-WFC 


6568/95 


300 


1 


1.3 


2010-02-09* 


2455237.37 


INT-fWFC 


6568/95 


1500 


5 


0.9 


2011-12-13* 


2455909.41 


INT+WFC 


6568/95 


2700 


9 


0.8 



Note. — The date indicates the start of observing night in the format yyyy-mm-dd. The date marked in boldface 
corresponds to the image that was used as reference for the astrometry and flux match of all images. Images from all 
dates marked with * were flux and PSF matched with the reference image. For 1999 only flux match was done due 
to its modest seeing. JD is the Julian Date at mid point of observations. In column 6 is the number of frames added 
together. 

^Published in lSeaquist et ahl i ll 9891 ) 
^Published in lShara et~an J2012D 
"^Published in lBode etldl J2004D 



Table 2 

Log of the spectroscopic observations 



Date 


Tclescope+ 


Slit P.A. 


Total exp. 


Number 


Resolution 


Spectral range 




Instrument 


(°) 


time (see) 


of frames 


(A) 


(A) 


2007-01-13 


INT+IDS 


22.5 


3600 


2 


1.3 


5730 - 7610 






48 


3600 


2 










315 


1800 


1 






2007-01-14 


INT-flDS 


69 


1800 


1 


1.3 


5730 - 7610 






293 


3600 


2 










334 


3600 


2 






2007-09-03 


NOT-fALFOSC 


31 


3600 


2 


0.7 


6350 - 6850 


2007-09-04 


NOT-fALFOSC 


86 


1800 


1 


0.7 


6350 - 6850 






107 


1500 


1 










312 


1800 


1 






2007-09-05 


NOT+ALFOSC 


49 


1800 


1 


0.7 


6350 - 6850 






173 


1800 


1 







Note. — P.A. is measured from North to East. Column 5 indicates a number of frames added together. 
See text for more details. 
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Some deviation from the circular symmetry is 
however visible because the diameter of the neb- 
ula in the NE-SW direction is smaller (105" as 
measured in the 2011 image) than in the NW- 
SE direction (118"). As described in details by 



Shara et al.l (|2012i ) , the nebula is composed of hun- 
dreds of knots and filaments of different sizes and 
brightnesses. Many of the knots have tails point- 
ing either toward or away from the central star. 
Shara et all (|2012[ ) counted 937 knots in their HST 
images. Some of them are blended in our lower- 
resolution ground-based images. On the other 
hand, the limited bandpass of the [Nil] F658N 
HST filter makes that some knots with the largest 
Doppler-shift displacements, which owing to the 
projection effects are located near the center of 
the remnant, are missed in the HST images. 

4. GEOMETRICAL DESIGNATIONS 

Throughout this article we use the designations 
as illustrated in Figure [51 In its left panel, we plot 
a view in the plane of the sky, where the line of 
sight corresponds to the Z axis and is perpendic- 
ular to the plane of the paper. The X axis points 
toward West, and the Y axis to North (up). Points 
1 and 2 indicate a knot's position at two different 
epochs. The displacement over one year (proper 
motion) is indicated by /i [arcsec yr~^], and its 
angle in the plane of the sky is a. PA indicates 
the position angle of the knot with respect to the 
central star as measured on our 20040129 frame. 
Both a and PA are measured from North to East 
counterclockwise. The apparent distance of the 
knot from the central star in the plane of the sky 
is named as d. 

When deprojecting the observed positions, each 
knot has a spatial distance i? from the central 
source (right panel of the Figure [2]) . The plane in- 
dicated in the figure is defined by the central star, 
the observer, and the knot considered: it is there- 
fore perpendicular to the XY plane, but changes 
depending on the knot under consideration. The 
speed of the knot is indicated by Ve2:p, and its com- 
ponents along the line of sight (radial velocity) and 
in the plane of the sky (tangential velocity) are in- 
dicated by Vrad and Vsfey, respectively. The angle 
between the line of sight and the velocity vector is 
d\. The angle that the line joining the central star 
with the knot forms with the line of sight is Q-i. in 




Fig. 2. — Designations used through this article. 
Left panel: sketch of a example knot as seen in 
the plane of the sky. Right panel: view in the 
perpendicular plane including the central star and 
the knot. See text for more details. 

the case of purely radial expansion, B\ = 02- Both 
angles are measured from the line of sight toward 
the velocity vector. 

5. PROPER MOTIONS 

Our multi-epoch imaging allows a precise deter- 
mination of the apparent expansion of the GK Per 
nebula in the plane of the sky. This is illustrated 
in the video in Figure |3] (only available in the elec- 
tronic version of the article). The animation in- 
cludes the early photographic and CCD images. 



• . Alt 

. ^"^2 . ... 
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Fig. 3. — Animation of the evolution of the 
GK Per optical nebula from 1953 to 2011. This 
animation is only available in the electronic ver- 
sion of the article. For arXiv version see ancillary 
file figure3_animation.mov . 
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Fig. 1. — Original frame from 20070905 obtained with the NOT. The image shows the filamentary nature 
of the GK Per remnant. The FOV is 2.8x2.8 arcmin^. 



astrometrically matched to the recent ones but in 
an arbitrary intensity scale and with a generally 
poor resolution, as well as a selection of the best 
INT images in the sequence of highly homogeneous 
data taken starting on 2004, in the same linear in- 
tensity scale. In the latter part of the sequence, 
not only the overall expansion of the outflow is 
clearly revealed, but also changes in brightness of 
some of the knots, as it will be discussed in the 
Section m 

The motions in the plane of the sky for 282 
individual knots were determined using the nine- 
teen, carefully registered images obtained between 



2004 and 2011. The measured knots are indicated 
in Figure The position of knots at each epoch 
was computed by fitting independently in the X 
and Y directions a Gaussian using the IRAF task 
imexamine, with fitting parameters adjusted to 
the mean characteristic of the knots. Only knots 
for which a reasonable Gaussian fit can be ob- 
tained were considered. Then the proper motion 
of each knot were determined by means of a least 
square fit of the positions at the different epochs. 
Both the distance variation from the central star, 
as well the separate motions in the X and Y direc- 
tions were computed, as illustrated in Figure[SJ In 



1990 1995 2000 2005 2010 




20" • 



Fig. 4. — Knots with proper motion measurements 
(crosses, blue in online version) indicated on our 
reference frame 20040129. Circles (blue in online 
version) show stars inside or nearby the remnant. 
FOV 2.8x2.8 arcmin^. 



this way, the proper motion vector is determined, 
which in addition to the magnitude of the spatial 
displacement /i also contains the information on 
its direction (angle a — aician{^x / l^-y) , see Fig- 
ure [2]). For all of the knots, a straight line pro- 
vides an excellent fit to the variation of d during 
the 7.9-yr lapse of time considered. This is further 
illustrated in Figure |6l Proper motions range from 
0".007yr-i to 0."53 yr^^. 

The prime focus camera of the INT has well 
defined and small geometrical distortion^ in the 
limited FOV covered by the GK Per nova rem- 
nant, and which are robustly removed during the 
astrometric registering of images. Therefore, er- 
rors in the proper motion determination mainly 
depend on the goodness of the fit of the position 
of the knots at the different epochs, which in turn 
depends on the knot's shape and brightness. As 
the proper motion error, we adopt the formal er- 
ror of the least square fit: its average value is 
0" .010 ±0" .006 yr-i. 

In the proper motion calculation, we did not use 
the images earlier than 2004 as they come from a 
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^http: / /www. ast.cam.ac.uk/~wfcsur/technical/astrometry/ 



-2000 2000 4000 6000 
Date [JD-2450000] 

Fig. 5. — Example of proper motion determination 
of a knot. The upper panel shows the variation of 
distance d with time, and the lower two panels its 
components along the X and Y axes {dx and dy) 
as defined in Figure [2j Open circles are data from 
1987 to 1999 (not used for the proper motion de- 
termination) , and filled circles those obtained from 
2004 to 2011. Dotted lines are the least square fits 
of data taken after 2004. Error bars are smaller or 
equal to sizes of symbols. 



less homogeneous set of observations with gener- 
ally poorer resolution or astrometric properties, 
as well as different filters. They are however used 
to search for signs of acceleration/deceleration of 
knots in the last decade with respect to the previ- 
ous 20 years. In this respect, we found no evidence 
for a systematic acceleration/deceleration, except 
in a few cases (see the lower three panels in Fig- 
ure [6]). Among these, a modest slow down is the 
most common case. 

The knots positions, their proper motions /z, 
/ia;, /ij^, and direction a, position angles, radial 
velocities, and the errors associated to all these 
quantities, are listed in Table [3l An excerpt is 



7 



5 


15 

10 

20 



15 - 



20 - 



15 ^ 



30 - 10.27 -26.09 



K 25h 
o 



50 - 



50 
45 
50 
45 
50 
45 



J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I |, 

4.11 0.38 



-12.24 0.56 



-8.34 13.19 



■ 21.8 -18.63 



a- 



10.68 -17.04 



;: 16.22 19.52 



..o.o- 



27.4 -24.92 



46.8 21.34 



: 3.72 -48.05 



-47.61 9.07 



30.06 38.02 



45 - 39.14 17.67 
40 

1 1 1 1 -I I'h ; I I I I 



o o 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I ll 

0.21 4.72 

o-« »•-►-»••**•■■■•- 



11.7 -0.16 



5 


H 15 

- 10 



■e-o- 



^m.'r-^. 30 

- 25 



...-o-o"- 



- 25 



,..-o'Cf ' 



20 
..■H40 

- 35 



-39.81 -12.3 



J: 52.2 -3.76 



:: 41.89 19.85 



19.32 30.49 



...-•■•o-o' 



— ■■■b' 



25.94 14.5 



0..Q- 



' ' ' ' 



45 
40 

55 
50 



- 50 



- 45 



40 

35 

30 
25 



-2000 



2000 4000 6000 -2000 2000 4000 6000 
Date [JD- 2450000] 



Fig. 6. — More examples of proper motion determinations. In the first seven rows, knots are ordered 
according to increasing proper motion magnitude (the range spanned by abscissae is the same for all graphs, 
10 arcsec). In the three bottom rows, some cases with signs of deceleration/acceleration in the last two 
decades are shown. At the upper-left corner of each panel, the knot's coordinates (dx and dy distances from 
the center, in arcsec, extrapolated at the reference epoch, see Table [3]) are indicated. Error bars are smaller 
or equal to sizes of symbols. See the caption of Figure [5] for the explanation of legends. 
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Table 3 

Knots positions, proper motions, and radial velocities 



(") 


dy^ 
(") 




(" s>— 1) 


(" »>— 1) 


(" v'^) 


(■■ V^^) 


(" yr-'^) 


PA 

(°) 


"PA 
(°) 


(°) 


(°) 


^rad 
(km 


0.06 


-27.43 


0.242 


0.005 


0.005 


0.004 


-0.242 


0.005 


180.1 


0.3 


181.2 


0.9 




0.09 


-48.91 


0.448 


0.006 


-0.009 


0.012 


-0.448 


0.006 


180.1 


0.2 


178.8 


1.5 




0.21 


4.72 


0.046 


0.011 


-0.007 


0.012 


0.046 


0.011 


357.5 


2.0 


8.9 


14.7 


808 


0.21 


4.72 


0.046 


0.011 


-0.007 


0.012 


0.046 


0.011 


357.5 


2.0 


8.9 


14.7 


919 


0.22 


33.44 


0.261 


0.003 


-0.005 


0.004 


0.261 


0.003 


359.6 


0.3 


1.0 


0.8 




2.25 


-9.25 


0.081 


0.004 


0.010 


0.010 


-0.080 


0.005 


193.7 


1.0 


186.8 


7.1 


-822 


3.41 


18.67 


0.153 


0.013 


0.003 


0.019 


0.155 


0.014 


349.7 


0.5 


358.9 


7.1 


721 


3.69 


-0.88 


0.048 


0.017 


0.041 


0.015 


-0.033 


0.022 


256.6 


2.5 


231.4 


21.0 


-894 


3.72 


-48.07 


0.497 


0.012 


0.057 


0.007 


-0.494 


0.012 


184.4 


0.2 


186.5 


0.9 




4.11 


0.38 


0.007 


0.006 


0.006 


0.006 


0.007 


0.011 


275.3 


2.3 


317.3 


54.0 


-794 



Note.— The entire tabic is publislicd only in the electronic edition of the article (For arXiv version sec ancillary file tabel3_online.txt). The first 10 lines 
arc shown here for guidance regarding its form and content. 

^Extrapolated to 2004-01-29, or ,1D2453034.37. 



presented in the paper, while the entire table is 
only available in the electronic version of the ar- 
ticle. In the table, individual knots are identified 
by their X and Y distance in arcsec from the cen- 
tral source {dx, dy) on the plane of the sky at the 
reference epoch as determined from the adopted 
least square fits. 

5.1. Comparison with published data 

A study of the proper motions i n GK Persei 



remna nt knots was recently done by iShara et al 



(|2012l ) using the 1995 and 1997 HST images listed 
in Table [H The proper motions of 938 knots is es- 
timated, but no quantitative c omparison w i th our 
measurements is possible as IShara et al. ( 2012 ) 
did not publish data for individual knots. In spite 
of the higher resolution of the HST, the seven- 
times longer baseline of our observations, the large 
number of measurements (19, which reduce statis- 
tical errors) , and the accurate knowledge of the as- 
trometric properties of the INT prime-focus cam- 
era, provide proper motion measurements which 
are mor e accurate by a f actor of four (cf. Fig- 



ure 8 in IShara et al.l |2012| with our average error 



of 0."01 yr-i). 

Anupama fc Prabhul (|l993h (AP93) have mea- 



sured the motion in the plane of the sky of 20 
individual bright knots using two ground based 
images obtained on 1984.63 and 1990.10 at Vainu 
Bappu Observatory (India), with a similar filter 
as in our observations. From visual inspection of 
their Figure 1, it is possible to safely identify 15 
knots in common with our measurements. The 
comparison of their proper motions is presented 



in Figure [T) Individual differences between the 
two datasets are above the quoted errors, but no 
systematic effects appear. Given that our mea- 
surements come from the analysis of 19 accurately 
registered high-quality images, we keep our proper 
motion measurements for the following analysis. 
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Fig. 7. — Compa rison of proper motion of 15 knots 
in common with lAnupama fc Pra ([1993). 



6. RADIAL VELOCITIES 

Our medium resolution spectra include the Ha 
line and the [Nil] doublet. No other emission lines 
were detected at the depth of these observations 
in the relatively small spectral range covered (see 
Tabled]). First, the spatial regions which clearly 
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correspond to knots identified in the images ob- 
tained at the same epoch were extracted from 
the long-slit data. In this way, 217 knots could 
be identified and measured. Then the signal in 
each region was integrated to create the Ha and 
[Nil] line profiles of each knot. They are ge nerally 



broad, as also noted bv IShara et al.l ()2012l ). with 



an instrument-corrected full width at half maxi- 
mum up to 200 km s~^, as illustrated in Figure [H 
The figure also shows that knots toward the center 
of the nebula (right panel) frequently have marked 
asymmetric shapes, with tails extending to veloc- 
ities as large as 300 km s^^. Given the strong 
projection effects in these central regions, it is not 
clear if these tails are intrinsic to the knots, associ- 
ated with related filaments, or caused by apparent 
superposition of more than one knot. 

The radial velocities of each knot were then 
measured by Gaussian fitting the [Nil] 6583 A 
line profile. When the line clearly had more than 
one peak or possessed non-negligible wings, multi- 
Gaussian fitting was applied, but the secondary 
components (which takes into account the wings 
or minor peaks) were not used in the following 
kinematical analysis. Errors on the radial velocity 
have been estimated by comparing multiple mea- 
surements of the same knots in different nights and 
with the two different telescopes used. A typical 
error of 8 km s~^ is adopted. All velocities were 
corrected to the Local Standard of Rest (LSR). 
As the LS R systemic velocit y of GK Per was esti- 
mated bv lBode et all (|2004l) to be 45 ± 4 km s-\ 
radial velocities were also corrected by this value. 
Table [3] lists the radial velocities of knots with re- 
spect to the systemic velocity of GK Per. They 
range from —989 to -1-967 km s^^. 

7. KINEMATICAL ANALYSIS 

Given its roundish overall morphology, as a first 
approximation we consider the nebula of GK Per 
to be a quasi-spherical shell in which knots expand 
radially from the central star. This hypothesis is 
broadly consistent with the plot of the observed 
radial velocities Vrad of the knots as a function of 
their apparent distance from the central star (Fig- 
ure m upper panel) . The two solid lines in the 
figure indicate the expected distance- velocity plot 
for spherical shells of radii of 35 and 55 arcsec, and 
expansion velocities of 400 and 1000 km s~^, re- 



spectively. Most of the knots are comprised within 
these two limits, suggesting that the knots of the 
nebula of GK Per form a relatively thick shell ex- 
panding with a significant range of velocities. 

This conclusion is supported by the behavior 
of the proper motions /i. The component of the 
expansion velocity in the plane of the sky imme- 
diately follows from /i according to equation (in 
convenient units) 



Vsky [km s' 



4.74-^ ["2/r-i]-i^ H, (1) 



where D is the distance to GK Per. For a distance 
of 400 ± 30 pc (Section [T]), Vgky of the measured 
knots varies between 13 and 1005 km s~^, and 
its variation with distance from the center is also 
consistent with a spherical thick shell model with 
parameters as described above (Figure ^ lower 
panel). 

The combination of the velocity components 
in the plane of the sky and along the line of 
sight provides the velocity vectors of the knots 
(see the right panel of the Figure [5]) once the 
distance to GK Per is fixed. We could measure 
both velocity components for 117 knots: their 

deprojected speed is Vexp = 



sky 



^racv and 



di = aiciaR{v sky /vrad) is the angle between the 
line of sight and the velocity vector (as in Fig- 
ure [5]). Assuming a purely radial expansion, 9i is 
also the angle that the line joining the central star 
with a knot forms with the line of sight, 62 in Fig- 
ure [2] The deprojected distance from the center 
would therefore be R = d/ \ sin^ij, and its projec- 
tion along the line of sight is Z = R cos 61. We 
define the sign of Z so that it is positive when a 
knot is approaching us (i.e. blue-shifted) and neg- 
ative when it is receding (red-shifted). We also 
define d to be positive for all knots in the North- 
ern side of the nebula and negative for those in the 
Southern half. With these notations, a represen- 
tation of the "depth" of the nebula along the line 
of sight (i.e. its extension along the Z direction) is 
shown in Figure 1101 In the figure, the observer is 
on the right side, and the plane Zd represents the 
folding of all planes perpendicular to the plane of 
the sky and passing through the central star and 
the knots considered (right panel of Figure [2]) . 99 
knots with a computed error on R smaller than 10 
arcsec are plotted. As stated above, the extension 
of the nebula along Z depends on the adopted dis- 
tance D. We assume that the distance value that 
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Fig. 8. — Example of intensity profiles of knots on the edges (left panel) and near the center (right panel) of 
the remnant. Vertical dotted line indicates the measured radial velocity for the corresponding knot. 



better corresponds to an overall spherical shape of 
the nebula is our "expansion-parallax" determina- 
tion of the distance to GK Per. The adopted value, 
estimated by visual inspection of the variation of 
the overall shape of the nebula as a function of D, 
is 400±30 pc. This is broadly equivalent to as- 
suming that Vsky for knots at the outer apparent 
edge of the nebula is the same as v^ad for knots in 
the innermost regions. 

Figure [Tl] shows the relations between depro- 
jected quantities, namely the speed Vexp and the 
distances R from the central star. As in previ- 
ous graphs, it confirms that the nova remnant is 
thick (nearly a half of its outer radius) and that 
the expansion speed of knots covers the signifi- 
cant range from 250 to 1100 km s~^, with the 
majority of them having a speed between 600 and 
1000 km s^^. In the figure, we also indicate the ex- 
pected relationship between velocity and distance 



for a ballistic ejection (Hubble-like flow) with an 
age t„oi;a=102.9 years, which is the time lapse 
between the nova maximum of 1901 Feb 22 and 
our reference epoch of 2004 Jan 29. Points above 
the solid line have a kinematical age smaller than 
tnova, and those below the line have a larger kine- 
matical age, which indicates that they suffered 
some deceleration in the course of the expansion. 
For comparison, in the figure we show as a dotted 
line also the relation expected for a ballistic ex- 
pansion with an age of 140 years. The knots age 
is further investigated in the next section. 

7.1. Kinematical ages 

A direct estimate of the kinematical age t of 
each knot can be obtained from the proper motion 
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distance d [arcsec] 

Fig. 9. — Above: Doppler-shift velocities of 217 
knots as a function of their distance from the cen- 
tral star projected in the plane of the sky. Solid 
lines indicate the position-velocity plot for thin 
spherical shells with the parameters indicated in 
the text. Below: the same, but for the velocity 
component in the plane of the sky from the mea- 
surements of the proper motions of 282 knots. 

measurements, according to 



t [yr] = 



(2) 



where d is the distance of the knot from the cen- 
tral star on the plane of the sky that we measured 
at our reference epoch (see Figure [2]). 

In Figure [T^ we present the kinematical age 
computed in this way as a function of the knots 
position angle. In the figure, only the 151 knots 
with an estimated error on the age smaller than 
8 years, and with an apparent distance from the 
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Fig. 10. — Nebular geometry along the line of sight 
{Z direction, observer on the right side) for the 
adopted distance to GK Per of 400 pc. Circles 
indicate inner and outer radii as in Figure [9] See 
text for details. 



centre d > 35 arcsec (in order to focus on the di- 
rections close to the plane of the sky and minimize 
projection effects) are included. Ages vary be- 
tween 96 and 170 years, the mean value, weighted 
by errors, being 118 ± 12 years. The horizontal 
line indicates the nova age tnova, and shows that 
in general knots have suffered only a modest de- 
celeration since their ejection. This conclusion is 
consistent with the comparison between the knots 
speeds of up to 1100 km s~^ determined above 
(Figure [TT|) , and an i nitial expansioii of the out- 
flow of --1300 km s"^ (|Pottaschlll959l ). Indeed, in 
the simplistic hypothesis of con stant deceleratio n 
throughout the knots lifetime ( Duerbec^ 1987 ). 
their observed mean age of 118 years implies that 
the present-day average speed is 77% of the orig- 
inal ejection velocity. The time ti for individual 
knots to decrease in speed by a factor of two would 
be 220 years. Even for knots with an apparent age 
of 140 years (the upper limit for the vast majority 
of points in Figures [TT] and [T2|) . velocity would 
have decreased by only 42% in the 102.9 years 
since the nova explosion, and ti would be 120 
years. Data therefore excl ude the short tim escale 
ti=58 years estimated bv iDuerbeckl (1987). 
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Fig. 11. — Expansion speed vs. linear distance 
from the central star for 13=400 pc. The dotted 
vertical lines indicate the radius limits of 35 and 
55 arcsec as in the previous graphs. The solid line 
indicates the relation expected for ballistic expan- 
sion of material ejected in 1901 (102.9 years from 
our reference epoch). The dotted line shows the 
same relation for an age of 140 years. 



The calculations of lDuerbeckl (|1987l ) were based 
on the measurement of the radius of the SW side 
of the nebula from images taken between 1917 
and 1984. The se measurements were extended by 
Seaauist et al. I ((19891 who analyzed them using a 
more physical model which considered the expan- 
sion of an isothermal or adiabatic shock created 
by the interaction of the expanding outflow with 
a tenuous circumstellar m edium. Such an analysis 
was fu rther extended by lAnupama fc Kantharia 
(|2005[ ) with images taken till 2003. We refrain 
from further extending the method, as we found 
that the measurement of the size of the nebula is 
highly subjective, owing to the difficulty in defin- 
ing the nebular edge due to the lack of an abrupt 
fall of surface brightness, dumpiness (which is en- 
hanced by the better resolution of the recent im- 
ages), and asymmetries. Indeed, our estimates of 
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Fig. 12. — Kinematical ages as a function of the 
knots position angle, for knots with d > 35". The 
horizontal line is the real age of the nova, 102.9 
years. 



the size differ from those of Seaauist et al. ( 19891 ) 
and Anupama fc Kantharial ( 20051 ) bv up to sev- 
eral arcseconds, which prevents any detailed dis- 
cussion or refinement of the results obtained by 
these authors. 

7.2. Deviations from spherical symmetry 

The proper motion measurements also allow us 
to measure any deviation from circular symmetry 
of the expansion pattern in the directions close 
to the plane of the sky. Anisotropics along the 
line of sight depends on the adopted distance of 
the system, and indeed in the previous section 
we used the argument inversely, having minimized 
any anisotropy in order to estimate a distance to 
the system. 

Considering the data in Figure [121 the mean 
weighted age in the different quadrants is: 125 ± 14 
yr (NE, 34 knots), 115 ± 11 yr (SE, 45 knots), 
113±10yr (SW, 34 knots), and 121±llyr (NW, 
38 knots). We conclude that no strong evidence 
exists for deviations from the circular symmetry 
in the expansion pattern. If any, knots in the 
Northern side of the nebula seem to have suffered 
a stronger deceleration throughout the remnant 
lifetime than those in the Southern part. There- 
fore, o ur calculations do not support the predic- 
tion of iBode et al.l (|2004l ) that knots in SW quar- 
ter should have been slowed down more than the 
rest. Smaller proper motions in the NW quadrant 
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Fig. 13. — Deviations of the proper motion veloc- 
ity vector a from the position angle PA of knots. 



are also mentioned in IShara et all (l2012h . 

In addition, our initial hypothesis of purely ra- 
dial motions can be tested by comparing the di- 
rection of the velocity vectors in the plane of the 
sky, a, with the position angle PA of knots in the 
plane of the sky, i.e. their radial direction from 
the central star. As in Figure [T^l only knots with 
d >35 arcsec are considered in order to limit pro- 
jection effects. Figure [T3l shows that velocity vec- 
tors are generally aligned along the radial direc- 
tion but a clear symmetric pattern of non-radial 
velocities is also observed. The maximum devi- 
ations are observed in the NE quadrant, where 
knots trajectories seems systematically bended to- 
ward the East with respect to the radial direction, 
and the Western part of the nebula, where knots 
velocity vectors are preferentially bended toward 
the South. 

8. BRIGHTNESS EVOLUTION 

The GK Per optical nebula suffered clear mor- 
phological and brightness changes in the last 50 
years. This is illustrated in Figure I14[ where a 
x2 zoom of the 1953 image is compared to a 
smoothed version of our 2011 image. The main 
changes are a relative brightening of the whole 
Eastern side of the nebula compared to the West- 
ern side, and in particular the development of a 
more prominent "bar" in th e NE edge (see also 
Anupama fc Kantharia 2005[) . and a progressive 
circularization of the SW edge of the nebula, i.e. 





Fig. 14. — Comparison of the 1953 photographic 
image of GK Per with a smoothed version of our 
2011 CCD image. The former is displayed with 
a x2 zoom with respect to the latter for easier 
comparison. 

of the region of the strong interaction observed at 
radio and X-ray wavelengths. In general, the out- 
flow looks more circular and uniform than in the 
past. 

As most of our images were taken with the 
same instrument, detector, and filter, the bright- 
ness variation of the nebula in the last decade can 
be quantified. The 15 INT-fWFC flux-matched 
frames marked with an asterisk in Table [T] were 
used. To measure the brightness of the nebula, all 
stars inside or nearby the remnant were masked. 
We first measured the brightness variation of the 
nebula as a whole. In Figure [TS] we present the 
evolution of the total fiux of the remnant normal- 
ized to the value that it had in 1999. It shows that 
the total iJa+[NII] fiux of the nebula has been lin- 
early decreasing, for a total of 31% over the 12.04 
years considered, i.e. at a rate of 2.6% per year, 
which is similar to what measured in the radio by 
Anupama fc Kantharial (|2005[ ) and attributed to 
adiabatic expansion of the remnant into the sur- 
roimding medium. 

Second, the fiux variation for each of the four 
quadrants (NE, NW, SE, SW) was computed. As 
can be seen from Figure [121 the NE quadrant 
shows the shallowest decline, with even some pos- 
sible re-brightening in the last 3 years. The quad- 
rant fading at the largest rate is the SW one. 

Third, we have also considered the fiux vari- 
ations for a number of individual knots. Their 
fiuxes were measured using circular apertures 
matched to the knots size, using an external an- 
nulus for background subtraction. The 1999 im- 
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Fig. 15.— Total i/a+[NII] flux of the remnant 
from 1999 to 2011. 
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Fig. 16. — Top panel: flux variations for different 
quadrants. Lower panel: flux ratio variations of 
NE and SW quadrants. 



age could not be used for the measurement of 
individual knots because of its poorer seeing. Af- 
ter excluding knots in crowded areas, we could 
measure the flux variation of 85 knots. Illustra- 
tive examples are presented in Figure 1171 In the 
figure, the flux is normalized to the maximum 
flux of each knot in the period considered. Out 
of 85 knots, 50% show a similar fading as the 
nebula as a whole (first two rows in the figure). 
Some 15% have a constant flux (third row), while 
some 20% of the knots show instead a monotonic 
brightening (fourth row in Figure [TTl) . Finally, 
15% of the knots are rapidly fading or brightening 
(or both) during the 2004-2011 period (last three 
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Fig. 17. — Illustrative example of the brightness 
changes of individual knots. The knots dx and dy 
coordinates as in Table |3] are indicated. 

rows). The most extreme cases show a sudden 
brightening of a factor of 5 in one year (bottom- 
left panel), or a similarly fast brightening and 
then a dramatic fading on a timescale of less than 
three years (bottom- right panel). The brightness 
variations o f indiv idual knots is also discussed by 
Shara et al.l (l2012h . 



9. DISCUSSION 
9.1. Distance 

Our kinematical determination of a distance of 
400±30 pc for GK Per (Section [71) compares well 
with the value derived using the absolute peak lu- 
minosity of novae. We have considered four dif- 
ferent relationships between My at p eak luminos- 
ity and the rate of the nova decline ( Cohen 1985; 
Capaccioli et al. 1989t della Valle fc Livio 1995 : 
Downes fc Duerbeckl 2000 ). depending on t2 (the 
time in days the nova takes to fade 2 mag from 
the peak brightness). Assuming the following for 
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GK Per: t2 = 7 d, a peak visual magnitude of 
m„= 0.2 maK, an d a reddening of E{B — V) = 0.3 



(|Wu et al.lll989h . we have computed the distance 
to GK Per using the different relationships. The 
average distance is 440±70 pc, in good agreement 
with our kinematic estimate. 

Our distance also broadly agrees with (and 
should improve) the estimates from other authors 
using the expansion parallax of the r emnant. The 
most w i dely used di s tance is 470 pc bv lMcLaughlin 



()l960l) . 

be 525 



Duerbecld (|1981r ) found the distance to 
pc. A m ore recent determination by 



Slavin et all (|1995^ states 455 ± 30 pc. 

Finally, a rather shor t distance to G K Per, 337 
pc, was determi ned in Warneii (|l987l) using the 
Bailey's method (|Bailevlll98l[) . This method con- 



siders the i^-band magnitude of the secondary 
star, which for cataclysmic variables is not eas- 
ily separated from the disc contribution: therefore 
the calculated distance should be considered as a 
lower limit. The largest distance de t ermin ation 
is the one bv lSherrington fc Jameson (1983), 726 
pc, using again the K mag and their estimated 
spectral type for the secondary star (K2 IV). 

9.2. Dynamics of the ejecta 

Our kinematical study, combining velocity mea- 
surements both in the plane of the sky and along 
the line of sights, has highlighted several ba- 
sic properties of the dynamical evolution of the 
GK Per nova remnant. The outflow consists of a 
thick shell composed of about a thousand knots 
and filaments. Knots velocities range from 250 to 
1100 km s^^, with the majority of them having a 
velocity between 600 and 1000 km s^^. The kine- 
matical ages of knots from the proper motion mea- 
surements indicate that most of them have been 
only modestly slowed down in the 102.9 years from 
the ejection to the present, with no significant de- 
celeration during the last 10-20 years. Clear ev- 
idence for systematic deviations of the expansion 
rate with direction is either not found. 

On the other hand, some long-term evolution 
of the overall characteristics of the nebula is high- 
lighted. The most notable changes are a progres- 
sive circularization of the SW nebular edge, and a 
relative brightening of the whole Eastern side com- 
pared to the Western part, with the development 
of a "bar" in the NE edge. 



In few words, the nebula seems to expand 
nearly ballistically, and is progressively becom- 
ing rounder and more uniform in terms of surface 
brightness. These results are somewhat surprising 
considering the previous interpretation of the evo- 
lution of the GK Per remnant. In fact, the asym- 
metry observed at radio and X-ray wavelengths 
(|Seaquist et al.lll989t iBalmanI lioosl ) seems to in- 
dicate that along the SW direction strong interac- 
tion of the remnant with relatively dense circum- 
stellar material has occurred. This is however not 
reflected in enhanced deceleration of the optical 
knots along this direction. 

The interpretation of the nebular dynamics 
should take into account these new results. A ba- 
sic role in the discussion is played by the origi- 
nal circumbinary density distribution into which 
the nova ejecta expanded. Its most visible expres- 
sion is the larg e bipolar ne b ula fir st detected from 



IRAS data bv iBode et al.l (|l987n . The strongest 



argument for the association between this bipolar 
structure and the nova is geometrical, that is, the 
emission is symmetrically placed around GK Per. 
The inherent problems of the classification of the 
bipolar nebulosity as an old planetary nebula from 
the p rimary star of the system were presented by 
Tweedy (1995};_ An alternative solution has been 



proposed bv lDoughertv et aL ( 1996f l: in it the sec- 
ondary star could be responsible for the formation 
of the extended nebula as it evolved off the main 
sequence triggering a major transfer of material 
into the primary star causing a born-again asymp- 
totic giant branch event. As revealed by deep opti- 
cal imaging the extended emission presents a typ- 
ical bipolar morphology with the polar axis ori- 
ented in the NW-SE direction. The nebula shows 
an asymme t ry bet ween the NE and SW edges that 



Bode et al.l (|2004l ) interpreted as the motion of 
the nebula to the SW assuming high relative ve- 
locity with respect to the I SM. Within this pic- 
ture the non-thermal radio [a nuDama & Prabh 
1993) and the X-rays (,Bal man & Ogclman 199 
emission in the SW direction is interpreted by 



i 



Bode et al.l (|2004^ to be a consequence of the nova 



ejecta ploughing the denser medium of the SW 
side of the bipolar nebula. 

We can test this idea in the framework of the 
existing hydrodynamical models of nebular inter- 
action with the ISM. The central binar y is moving 
with a space velocity of 45 ± 4 km (jBode et al 
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20041 ) ■ or 4.6xl0~^ pc yr~^. For simplification, 
we consider the central sourc e as a sing le star (as 
successfully done for Mira, see lMartin et al. 2007). 
Moreover, given that the direction of the move- 
ment is almost perpendicular to the axis of the 
bipolar lobes we can also assume spherical symme- 
try (the formation of a bipolar structure is com- 
plex and beyond the scope of the present work). 
With this simplifications at hand we can consider 
that the interaction process can be similar, at least 
regarding the discussion that follows, to that of 
the wind of an AGB star with the ISM (see e.g. 
Villaver et al. I l200,l I2OI2I ). In particular a simu- 
lation wi th a relative speed o f 50 km s~^ can be 
found in I Villaver et al.l (|2012l) . These simulations 
show that once a bow-shock structure is formed, 
the subsequent stellar wind expands unperturbed 
inside the asymmetric outer cavity created by the 
interaction. In these conditions, never the stel- 
lar material ejected inside this cavity can feel the 
outer interaction if it is inside the bow-shock struc- 
ture. The same idea could be applied to GK Per: if 
both the outer bipolar nebula and the nova ejecta 
are the result of the evolution of the central binary, 
and the nova shell is expanding inside the bipolar 
nebula, then the ejecta cannot see the interaction 
with the ISM, and should expand unperturbed. 
This interpretation is consistent with the lack of 
evidence for deviations of the expansion rate with 
direction found in this work. In particular, the 
SW side, where most of the interaction with the 
circumstellar medium was assumed to occurred in 
the past (see radio and X-ray data) does not show 
any sign of larger deceleration than the rest of the 
nebula. The asymmetry observed in the radio and 
X-rays, however, would remain unexplained in this 
scenario. 

An alternative possibility is the existence of a 
mechanism that allows interaction with the ISM 
to permeate from the outer structure to the nova 
shell. An interaction wi th magnet ized ISM could 
provide this mechanism |Sokei;_&_Dgani 1997). In 
fact, the nature of the radio and X-ray emission 
could be explained as dense gas clouds moving 
supersonicall y th rough a magnetized low-density 
medium. I Jones et al.l (jl993 ) show how the syn- 
chrotron emission from relativistic electrons in- 
creases rapidly when the cloud begins to fragment 
under a Ray leigh- Taylor instability. They also 
show how the interaction of this wind with dense 



condensations can lead to strong magnetic fields. 
If relativistic particles are formed as the fast wind 
is shocked, then the enhanced magnetic field will 
result in non-thermal radio emission. 

At this point we can turn our attention back 
to the outburst mechanism that leads to the for- 
mation of the nova ejecta. Its thick shell geome- 
try can be the result of an ejection with a range 
of initial velocities, or the dynamical evolution 
of the interaction of different components. We 
favour the latter hypothesis, which is supported 
by the observations of nova shells soon after out- 
burst. It has been argued in the literature that 
two distinct outfiow components of different ori- 
gin are present at early times: an earlier, slow one 
from the circumbinary gas, and a later, fast one 



result of the outburst (see e.g . ICassatella et al 
l2004t IWilhams fc MasonI I2OI0I) In t he case of 
Nova Cvg 1992. ICassatella etM (|2004[ ) show that 
the main shell ejected at the moment of the out- 
burst contains the bulk of the ejected matter and 
is located in the outermost region of the ejecta re- 
gion, while a less massive impulsive event forms a 
faster expanding shell that eventually will catch up 
with the main shell generating a region of shock 
interaction where the two shells meet. Eventu- 
ally, the process gives rise to the formation of a 
new shell formed by the merger of the two. A 
similar scenario in terms of d iscrete sh ell ejection 
events is proposed bv Wil lia ms et aP (2008) and 
Williams fc MasonI (|2010l ) to explain the absorp- 



tion line spectra of the early ejecta of 15 nova 
shells. 

In GK Per we are probably witnessing the late 
evolution of such of a process, and in this respect it 
is natural that any asymmetry originally present 
in the dense eje cta (see e.g. the 1917 images in 



Bode et al.ll2004l ) disappears with time as the com- 
bined shell expands outward. This explains the 
progressive circularization and homogenization of 
the optical outfiow. The original asymmetry can 
also be understood in te rms of the binary i nterac- 
tions in the system (e.g. Svtov et al. 20091) . 



Based on observations made with the Isaac 
Newton Telescope, operated by the Isaac New- 
ton Group, and with the Nordic Optical Tele- 
scope, operated jointly by Denmark, Finland, Ice- 
land, Norway, and Sweden, on the island of La 
Palma in the Spanish Observatorio del Roque 
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de los Muchachos of the Institute de Astroffsica 
de Canarias. INT data were in part reduced 
by the lAC support astronomer staff. This pa- 
per also makes use of data obtained from the 
Isaac Newton Group Archive which is maintained 
as part of the CASU Astronomical Data Cen- 
tre at the Institute of Astronomy, Cambridge. 
This work makes use of EURO-VO software, Al- 
adin and TOPCAT. The EURO-VO has been 
funded by the European Commission through con- 
tracts RI031675 (DCA) and 011892 (VO-TECH) 
under the 6th Framework Programme and con- 
tracts 212104 (AIDA) and 261541 (VO-ICE) un- 
der the 7th Framework Programme. Data pre- 
sented here is partially based on photographic 
data of the National Geographic Society - Palo- 
mar Observatory Sky Survey (NGS-POSS) ob- 
tained using the Oschin Telescope on Palomar 
Mountain. The NGS-POSS was funded by a grant 
from the National Geographic Society to the Cal- 
ifornia Institute of Technology. The plates were 
processed into the present compressed digital form 
with their permission. The Digitized Sky Sur- 
vey was produced at the Space Telescope Science 
Institute under US Government grant NAG W- 
2166. The compressed files of the "Palomar Ob- 
servatory - Space Telescope Science Institute Dig- 
ital Sky Survey" of the northern sky, based on 
scans of the Second Palomar Sky Survey are copy- 
right (c) 1993-2000 by the California Institute of 
Technology and are distributed herein by agree- 
ment. This research has made use of the USNOFS 
Image and Catalogue Archive operated by the 
United States Naval Observatory, Flagstaff Sta- 
tion (http: / / www. nofs. navy, mil/dat a/fchp ix/ ) . 
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